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ABSTRACT 



Aims. We searched for a fast moving Htt shell around the Crab nebula. Such a shell could account for this supernova remnant's 
missing mass, and carry enough kinetic energy to make SN 1054 a normal Type II event. 

Methods. Deep Ha images were obtained with WFI at the 2.2m MPG/ESO telescope and with MOSCA at the 2.56m NOT. The data 
are compared with theoretical expectations derived from shell models with ballistic gas motion, constant temperature, constant degree 
of ionisation and a power law for the density profile. 

Results. We reach a surface brightness limit of 5 x 10"* erg s^'cm^^sr"'. A halo is detected, but at a much higher surface brightness 
than our models of recombination emission and dust scattering predict. Only coUisional excitation of Ly/3 with partial de-excitation 
to Ha could explain such amplitudes. We show that the halo seen is due to PSF scattering and thus not related to a real shell. We also 
investigated the feasibility of a spectroscopic detection of high- velocity Ho- gas towards the centre of the Crab nebula. Modelling of 
the emission spectra shows that such gas easily evades detection in the complex spectral environment of the Ha-line. 
Conclusions. PSF scattering significantly contaminates our data, preventing a detection of the predicted fast shell. A real halo with 
observed peak flux of about 2 x 10"^ erg s 'cm^^sr^' could still be accomodated within our error bars, but our models predict a factor 
4 lower surface brightness. 8m class telescopes could detect such fluxes unambiguously, provided that a sufficiently accurate PSF 
model is available. Finally, we note that PSF scattering also affects other research areas where faint haloes are searched for around 
bright and extended targets. 

Key words. ISM: Supernova remnant, Supernovae:individual(SN1054) 



1. Introduction 

According to historical records, SN 1054 was visible dur- 
ing daytime for more than thr ee weeks and for almost two 
years during night-time (see e.g. Clar k & Stephenson|[T977l) . Its 
lightcurve (e.g. [Soller man et al. 200l|), hydrogen-rich nebula, 
and the presence of a centra l neutron star (IStaelin & ReifensteinI 
Il968t IComella et al.l Il969h are all compatible with a normal 
core-collapse event. However, the velocities o f the filament s in 
the Crab nebula are rather low (~ 1400 km s"','Woltier"1972), as 
is the observed ejecta mass of 4.6 +1.8 Mq (Fesen et al. 1993). 

The elemental abundances and the low velocity and mass of 
the filaments may be accomplished by a low-energy explosion 
of a 8 - 13 Mo star (Nomoto 1985, 1987; Kitauraet al. 200^. 
Indeed, supernovae with very low kinetic e nergies are known , 
such as SN 1997D with (1 - 4) x 10^" erg (iTuratto et al.|[T998t 



IChugai & Utrobinll2000l) . However, no events with as little ki- 
netic energy as SN 1054 (1 .5 x lO'*^ erg) have yet been identified. 

IChevalieil (119771 1 1985h suggested that a hydrogen-rich shell 
of a few solar masses could carry most of the kinetic energy. 
Such a shell could be observed as extended optical emission, or 
as non-thermal radio and thermal X-ray emission caused by a 
shock front advancing through the inter stellar medium such as 
in the case of SNR G2 1.5-0.9 (in X-rav. iMatheson & Safi-Harbl 
I2OO5I) . 

Searches for such emission aroun d the Crab nebula have 
been performed at rad i o wavelengths jWi lson & Weiler 19821; 
Velusamvl [19841 [1981 iTrushkinI [l986t [Velusamv et al. .19921 



samvl I 

Frail etal.1 1 1995b in t he optical (e.g iGuU & FesenI 119821 
Fesen & Ketelsen' '1985) and in X-rays (Mauche & Gor ensteinI 
,1985, 1989; Predehl & Schmitt 1995; Seward et al. 2006), yield- 
ing negative results. Possibly the Crab nebula expands into a lo- 
cal bubble of the ISM wi th correspondingly low interaction rates 
(see [Wallace et al.ll 19991 and references therein). 

The first ev i dence for the existence of a shell was reported 
by lHester et al.l (Il996l) . in the form of [O III] emission surrou nd- 
ing the outer boundaries of the nebula. ISankrit et al.l (1 19981) in- 
terpreted this emission as due to radiative shocks at the outer 
edge where the synchrotron n ebula accelerates the hypothetical 
shell. [Sollerman etlTI ( [20001) detected a blueshifted absorption 
feature in C IV /11550 in a far-UV HST spectrum, extending 
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Table 1. Summary of the WFI, MOSCA and FORSl observa- 
tions 



Instrument 


Exp. time [sec] 


Seeing 


Observation 




(number of exp.) 


["] 


dates 


MOSCA 


20700 (23) 


I'.'O 


2006-12-20 


WFI 


23520 (36) 


I'.'l 


2000-11-06/07 








2002-02-02/05/06 



up to ~ 2500 km s~'. Suc h a feature of the shell was predicted 
bv lLundqvist et al] (Il986h using photoionisation models. Lower 
limits of 0.3 Mq and 1.5 x 10"*^ erg were derived for the mass 
of the shell and its kinetic energy, respectively. This is signifi- 
cantly below the canon ical value of 10^' erg, but the results of 
ISoUerman etlTI (|2000|) do not rule out a normal Type II explo- 
sion. For a comprehensive summary of the search for a halo and 
for a general overview of the Crab nebula see ,Heste7](l2008l) . 

By means of sufficiently deep Ha surface brightness maps 
strong cons traints can be put on the shell mass and its radial den- 
sity profile. [Sollerman et al.l(l2000l) showed that a fast shell with 
4 Mq would escape the deepest Ha search so far, rea ching a sur- 
face b rightness limit of 1.5 x 10"^ erg s^'cm^^sr"' dFesen et alj 
[19971 using spectroscopy). In this paper we describe our at- 
tempts to image an Ha halo with the Wide Field Imager 
(WFI) (iBaade et alj|l999l) at the 2.2m MPG/ESO telescope and 
the Mosaic Camera (MOSCA) at the 2.56m Nordic Optical 
Telescope, reaching a sensitivity of 5.0 x 10"** erg s^'cm^^sr"'. 
The observations and data reduction are summarised in Sect. 2. 
We analyse the images in Sect. 3 and link the results to theoret- 
ical models in Sect. 4. The feasibility of a direct spectroscopic 
detection towards the centre of the Crab nebula is evaluated in 
Sect. 5. Our results are discussed in Sect. 6, followed by our con- 
clusions. Throughout this paper we use the term shell for actual 
physical material around the Crab nebula, and the term halo for 
its appearance in images and spectra. 

2. Observations and data reduction 

2.1. Imaging data 

WFI covers 34' x 33' with a pixel scale of 0'.'238. We observed 
in the Ha filter (A^ = 6588 A, FWHM=74 A) using a very 
wide dither pattern (see Fig. [1]) in order to suppress instrumen- 
tal effects and to obtain a good superflat. This is especially im- 
portant when sear ching for low surface brightness features (see 
lErben et aLll2005L and in particular their Fig. 8). Observations 
were carried out in service mode in dark and grey time (see Sect. 
IZTH for details). 

Independent Ha images (^ = 6564 A, FWHM=33 A, NOT 
filter code #21) were taken with MOSCA in dark time. MOSCA 
has 4 CCDs and offers in this filter a comparatively small, 5'.5 
diameter circular field with 0'.'217 per pixel. It was centred on 
the north-western part of the Crab nebula (see Fig.|2l). 

2.2. Data Reduction 

We proc essed the multi-ch ip camera data with the THELf] 
pipeline dErben et"an l2005l) . The images were overscan cor- 
rected (not for MOSCA, which does not have an overscan), debi- 
ased, flatfielded, gain-corrected, and in case of WFI also super- 
flatted. The MOSCA data could not be superflatted as the target 

' Available at |http://www.astro.uni-bonn.de/~mischa/theli.htniI| 
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Fig. 1. The WFI detector layout and the dither pattern. The Crab 
nebula was placed on three main positions on the detector and is 
represented to scale by the solid circles. At each of these main 
positions random dither offsets of up to +2'.5 were applied (small 
dashed circles). 



covered most of the field of view. One of the MOSCA CCDs 
(#4) did not flatfield properly and was partially masked (see also 
Fig. [3]). The images were then individually weighted and we 
determined the relative photometric zeropoints between them. 
After sky subtraction and a full astrometric calibration the im- 
ages were coadded. The photometric zeropoint is constant across 
the stacked mosaic images with an error of ~0™07. 

In the following we describe the superflatting and sky sub- 
traction. These two steps are critical as they must not suppress 
or enhance a halo while correcting for flat field residuals and 
inhomogeneous sky background. 

2.2.1. Superflatting 

Three main dither positions (Fig. [U were chosen for the WFI 
observations. Dither offsets of 4' - 5' were made at each one. 
Since the Crab nebula is 6' in size, this approach ensured that 
every pixel saw the sky much more often than the object, which 
is the basic requirement when mak ing superflats. We first run 
SExtractor jBertin & Arnoutslll996l) to detect all objects con- 
sisting of at least 5 connected pixels, each of which having a 
minimum S/N-ratio of 2. In this manner stars and the Crab neb- 
ula were masked out. Lower detection thresholds would have 
removed sky features as well which must be preserved. The 
masked images of a given CCD were then scaled to the same 
mode and median-combined using a 2.5cr outlier rejection. This 
is the standard approach in THELI. 

1 1 of the 36 images were taken in grey conditions, with the 
moon 12 degrees or less above the horizon and 80 to 110 degrees 
away from the target. 7 of these exposures were significantly 
affected by uneven sky background or reflections and removed 
from further processing. In the remaining 4 grey exposures only 
one or two of the chips were affected (none with the Crab nebula 
inside) and not used for the superflat. 
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The superflats were then smoothed with a 1 ' wide kernel in 
order to suppress small scale noise. Since CCDs 1 and 8 still 
showed significant, time-variable background features in all ex- 
posures, they were excluded at this stage. 

2.2.2. Sky subtraction 

After superflatting and rejecting CCDs 1 and 8, only small 
(1% - 3%) and large-scale (~ 10' ) gradients were left in the 
exposures. They are larger in size than the Crab nebula and a 
possible halo. The outermost 2' of the right edges of CCDs 4 
and 5 are occasionally brighter by about 2%. This is not a cause 
of concern as the dithering never moved the target close to this 
region. 

THELI uses a similar approach for background subtraction 
as for the superflat. All objects with at least 5 connected pixels 
and 2cr above the sky noise were replaced with the mode of the 
remaining pixels. This intermediate image was then smoothed 
and subtracted. We tested two different smoothing scales of 1'.8 
and 3'.0, and as an additional check also created a coadded im- 
age where only a constant sky was subtracted from each CCD. A 
very similar halo around the Crab nebula is detected in all three 
coadded images, thus it was not artificially introduced by the 
data reduction. In fact, we show below that it is caused purely by 
PSF scattering. The coadded images with the r.8 and 3'.0 back- 
ground smoothing scales are very similar, and we conservatively 
chose the latter for further analysis. A comparison with the con- 
stant sky subtraction image shows that the modelling lowered 
the halo's brightness by less than ~ 10%. 

The sky background of the MOSCA data was not modelled, 
as only small areas of sky were available between the Crab neb- 
ula and the edge of the field. We subtracted constant estimates 
obtained from empty areas located ~ 2' from the nebula's outer 
edge. Some particular issues with the MOSCA reductions are 
explained in the Appendix. 

2.2.3. Absolute photometric calibration 

We used the Crab pul sar and its flux-cahbrated spectrum (from 
ISollerman et alj|2000h for the absolute photometric calibration 
of the WFI data. Since the pulsar lies within a complex nebular 
environment, we must remove the contribution of the latter from 
the pulsar's flux. We obtained a DAOPHOT PSF model from 
isolated, bright but unsaturated, field stars and subtracted it from 
the pulsar's PSF. In this way the nebular contribution could be 
quantified and removed. The pulsar's magnitude was then aper- 
ture corrected, using the same field stars. The differences be- 
tween aperture and PSF photometry were then < 0.01 mag, and 
the final photometric zeropoint became ZPwfi = 21.42 + 0.15. 

Flux-calibration for the MOSCA data was based on 
DAOPHOT aperture photometry of the spectro-photometric 
standard star HD 93521. After appHcation of the MOSCA filter 
transmission curve to the reference spectrum, the photometric 
zeropoint became ZPmosca = 21.5 + 0.10. With this calibra- 
tion there was a constant offset of 0.07 mag in field stars seen by 
MOSCA as compared to WFI. The Carlsberg Meridian ArchiveQ 
shows increased extinction levels for the corresponding week in 
La Palma, however without a measurement for our night. The 
observations were thus probably not made in photometric con- 
ditions, and we corrected the MOSCA zeropoint for this offset. 



^ http://www.ast.cam.ac.uk/~dwe/SRF/camc_extinction.htmI 




Fig. 2. Part of the coadded WFI image. For clarity we did not 
invert the centre of the nebula. A faint halo around the SNR is 
visible. Superimposed is the slit position for the FORSl spectra, 
together with the location where we tested the feasibilty of a 
spectroscopic detection of a fast shell. The large circle shows 
the MOSCA field of view. 




Fig. 3. The coadded MOSCA image. The missing part to the up- 
per right is due to a partially masked CCD. 



3. Observational results 

In this section we present the results from our imaging cam- 
paigns. We trace clear but inconsistent signals in both the WFI 
and MOSCA data and discuss them in Sect. |4] 

Our halo detection strategy is described in the following 
based on the WFI data. It was done in a similar way for MOSCA. 
We 
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- mask all stellar flux so that it does not contribute to the halo, 

- define an outer boundary of the Crab nebula, 

- measure and average the halo brightness as a function of dis- 
tance from this boundary in many different sectors. 

3. 1 . Masking the image 

In order to mask any unwanted light, we ran SExtractor with a 
detection threshold of 2cr and a minimum number of 15 con- 
nected pixels. The noise in coadded mosaic images is very un- 
even due to the dither pattern and the gaps between the CCDs, 
and this fact must be taken into account in the detection process. 
The coadded weight map reflects these noise properties and was 
therefore slotted into SExtractor as well. In this way stars, filter 
ghosts, and the Crab Nebula itself were masked while the halo 
was preserved. Some residuals of scattered light were still visi- 
ble around brighter stars and masked manually. These masks are 
shown as black regions in Fig. |5] Also shown in this Figure is 
the contribution of PSF scattering mimicking a halo, which we 
discuss in Sect. 14. 11 

3.2. Defining tiie boundary around tlie Crab nebula 

We measure the halo brightness as a function of distance from 
the nebula in several sectors protruding outwards from the ex- 
plosion centre. Once the sectors cross the edge of the nebula, the 
radial coordinate (distance) is set to zero, and we begin integrat- 
ing the halo flux. However, the edge's fractal-like appearance 
(Fig. IHi makes the determination of the zero-point of the radial 
coordinate difficult. Isolated or half-connected groups of pixels 
still reside inside, leading to an underestimation of the distance 
between the explosion centre and the outer edge. As a result, the 
stack of individual profiles is not properly aligned, leading to in- 
creased scatter in the result. To avoid this, we defined the manual 
mask indicated by the solid line in Fig.|4] 

3.3. Tiie radial halo profile 

We assume the centre of the halo to coincide with the ex- 
plosion centre, located at Q'y2000 = 05''34"'32.7'', 5y2000 - 
-h2 2°00'48"7 (back-pr ojecting the pulsar proper motion as given 
by iKaplan et al.ll2008[) . The exact position is not a critical pa- 
rameter for our analysis, the measured halo profile changes neg- 
ligably when the assumed explosion centre is shifted by a few 
arcseconds. We defined 72 sectors, each 5 degrees wide, emerg- 
ing from this centre. Those sect ors that covered the jet-like Crab 
chimney (see lRudie et al.ll2b08L for details) were excluded from 
our analysis, hence 68 (MOSCA: 20) sectors remain. 

In order to scan the profile in each sector, we divided them 
into 29 (MOSCA: 19) equidistant radial bins 24 pixels in length 
corresponding to 5'.'7 (5'.'2) for WFI (MOSCA). Their mean 
brightness was calculated using an interative 2.5cr outlier rejec- 
tion. This clipping procedure reduced the error bars of the final 
radial profile by just 5%, indicating that our masks are already 
very efficient. 

Figure|6]shows the results. The error bars represent the stan- 
dard deviation of the mean surface brightness and are domi- 
nated by brightness fluctuations as a function of azimuthal an- 
gle. We can trace the haloes to distances of V9 and 13 for 
WFI and MOSCA, respectively. The MOSCA signal peaks at 
2.5 X 10"^ erg s"'cm"-sr"' and is twice as high as that for WFI, 
but it also decays faster. If x is the distance from the outer edge 
of the nebula in arcminutes, then the profiles are proportional 




Fig. 4. Defining the boundary (mask) of the Crab nebula. Only 
the south-eastern area is shown. The automatic masking with 
SExtractor left a fractal-like border which hampers the determi- 
nation of the starting point for the halo measurement in a sec- 
tor. The solid line shows the manual re-definition of the mask. 
Most of the stellar flux was automatically masked, but remain- 
ing residuals around brighter stars had to be covered manually. 
These masks can be seen in Fig.|5] 

to lo-o-^'*' (WFI) and lO-' '^^ (MOSCA). These inconsisten- 
cies are investigated in Sect. H) The sensitivity limit reached is 
5 X 10"*^ erg s^'cm^^sr"', corresponding to 31.2 mag per square 
arcsecond. 

4. Analysis 

We show here that the haloes are caused by PSF scattering, and 
we infer upper limits to the surface brightness of a real shell. We 
present model calculations for three different physical processes 
that can lead to Ha emission in the shell and compare them to 
our observations. 

4.1. PSF scattering 

Figure |7] shows filter ghosts and other scattered light around 
stars. In order to evaluate how much this effect contributes to 
the haloes, we obtained azimuthally averaged PSF models for 
both data sets (see Fig.O. 

The PSF model for WFI extends to a radius of 2'.7 and was 
created from very bright stars in the coadded image. The satu- 
rated core was replaced by one obtained from unsaturated stars. 
During the construction of the model a strict Icr rejection thresh- 
old was applied to suppress contamination from field stars. The 
resulting PSF profile exhibits several bumps due to the filter 
ghosts which make simple mathematical fits inadequate. We also 
created a second PSF model with reduced wing amplitude, by 
subtracting the RMS value of the sky background on a large 
scale from the first model. In this way we could evaluate how 
sensitive light scattering is to small changes in the amplitude of 
the PSF wings. 
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Fig. 5. Masks for the field stars, together with two of the 72 mea- 
surement sectors. These converge towards the explosion centre 
and are drawn to scale. The density of the radial bins is four 
times higher than indicated, and the field of view is similar to 
that in Fig. [J] The brightening in the immediate surrounding of 
the Crab nebula shows the amount of PSF scattering expected 
forWFI(see Sect.gl]). 




.0 1.5 2.0 

r[arcmin] 

Fig. 6. The halo profiles observed. They fall off following power 
laws with slopes of -0.88 (WFI) and -1.12 (MOSCA). 



The PSF model for MOSCA was obtained in a similar way 
from two moderately bright stars in a single 600s blank field 
exposure taken in the same night. The model's empirical part 
extends out to 0'.85, beyond which it is extrapolated by a power 
law with index [3 - -2.4 + 0.4. The slope was determined from 




Fig. 7. Illustrates the WFI filter ghosts. The star shown is the 
brightest seen in the lower left edge of Fig. |2] The core of the 
PSF is labelled as 1, followed by three filter ghosts of decreasing 
intensity, labelled 2-4. 



the profile between 0'.45 and 0'.85. Note that the uncertainty of 
this model at a radius of 0'.2 is already larger than that for the 
WFI PSF at 2'.0. 

To evaluate the PSF contribution, we zeroed all pixels in the 
coadded images having values less than 5 times the sky noise, 
keeping only the cores of stars and the Crab nebula. These im- 
ages were then convolved with the PSF models, and the profile 
measurements were repeated with identical settings and masks. 
The result for WFI is shown in Fig. |9] (dotted line), together 
with the contribution from the PSF model with lowered wing 
amplitude (dashed line). The latter significantly underestimates 
the halo flux as compared to the unaltered model. Residual light 
from field stars contributes to the halo only for radii larger than 
I'.5. We evaluated this by convolving an image where all field 
stars were masked out. Within the error bars of the best-fitting 
PSF and of the measurement, a real halo with an Ha peak flux 
of 2 X 10"^ erg s^'cm^^sr"' can still be accomodated, compara- 

erg s^'cm^^sr"' obtained by 



iFesen et" 



; upper lim 



Figure [TOl shows the corresponding result for MOSCA. The 
halo is very well explained by PSF scattering (dotted line). The 
other two lines represent the Icr uncertainty of the power-law fit 
to the PSF wing. Due to the less reliable MOSCA PSF model we 
could not improve the upper flux limit obtained with WFI. 

We stress that simplified approaches such as a comparison 
of the slope of the PSF with the slope of the halo, or of the 
brightness ratios between inner and outer part of the PSF and 
inner and outer part of the Crab nebula, are insufficient. In these 
ways PSF scattering would have been underestimated by about 
3 orders of magnitude. Only full 2D convolutions can quantify 
the effect accurately. It appears plausibl e that an earl ier detec- 
tion made in H a on a photographic pl ate (Murdin & C lark 19811; 
lMurdin|[T994b also shows this effect. [Davidson & Fesea (,19851) 
argued that the detection made was of instrumental nature, but 
the subject was not investigated further 

4.2. Theoretical models for Ha emission from a fast sfiell 

Our data still allows for a peak Ha surface brightness of 2 x 
10"^ erg s 'cm^^sr ' through the WFI Ha-filter . The u pper limit 
of 1.5 X 10"^ erg s^'cm^^sr"' by iFesen et al ] (Il997h refers to 
the total unfiltered surface brightness. We compare these val- 
ues below with predictions based upon the shell models of 
Sollermanet al. (2000). Therein, the shell is of spherical shape 
and expands freely outside the Crab nebula with V cc R. It has 
an inner minimum radius of R„iin = 5.0 x 10"^ cm corresponding 
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Fig. 8. The two PSF models for WFI (soHd line) and MOSCA 
(dotted line). The bumps in the WFI PSF at 0.15, 0.4 and 0.8 
arcmin reflect the various filter ghosts. The MOSCA PSF is ex- 
trapolated for radii larger than OISS by a power law with index 
jS - -2.4 (not shown). MOSCA has a significantly more com- 
pact PSF than WFI for radii larger than about 4" (0'.06), which 
explains the steeper slope of the halo in that data. 
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Fig. 9. The WFI halo profile with the contributions from the PSF 
(dotted line), and from the PSF with reduced wing amplitude 
(dashed). 



to a minimum velocity of 1680 km s"'. Unlike the Crab nebula 
the shell is not affected by the pulsar wind. The mass of the shell 
is set to Msh - 4.5 Mq, it has a fixed electron temperature of 
Te —2x 10"* K dLundqvist et al.| [l986), and a radial density pro- 
file p oc R '^, with /; = 3 or 7/ = 4. Sma Uer values for rj do no t 
fit the profile of the C IV absorption (see lSollerman et al]|2000h . 



and values rj > 5 would make the Ha surface brightness surpass 
our detection limit. 

All fluxes determined from the models must be corrected 
for galactic extinction before they can be compared to obser- 
vations, using the reddening pa rameters E{B - V) = 0.52 and 
R = 3.1 dSoUerman et al.ll2000l) . Our dereddened upper limit of 
the Ha surface brightness through the WFI Ha-filter becomes 
6.6 X 10"^ erg s"'cm"-sr"'. 



4.2.1. Recombination emission 

For the modelling of the recombination emission we furthermore 
assume that hydrogen in the shell is fully ionised and that the gas 
has a He/H-ratio of 0. 1 by number Figure[TT]shows the modelled 
surface brightness profile as a function of distance from the edge 
of the nebula for two cases: rj = 3 (top panel) and // = 4 (bottom). 
The solid lines describe the surface brighness captured by the 
WFI Ha-filter. 

For comparison we also plot the signals expected through 
the WFI red continuum filter (dashed line) and for unfiltered 
observations (dotted line). The continuum filter selects Ha- 
emission from shell gas receding from the observer with 2000 - 
6000 km s"'. The advantage is that the Crab nebula is much 
less bright in these wavelengths and thus the amount of scat- 
tered light is reduced, which favours a detection. On the other 
hand, the gas seen in this filter is actually at significantly larger 
distance behind the Crab nebula and just projected next to the 
nebula's edge. The densities in these parts of the shell are much 
smaller and thus the emission is significantly reduced (see Fig. 
fTTT i. The dotted line represents the fluxes expected for unfiltered 
observations which we need later on for comparison with spec- 
troscopy. 



Peak surface brightness for the Ha-filter occurs at the 
edge of the nebula with 1.5 x 10"^ erg s"'cm"-sr"' for rj — 
4. For comparison, the observed WFI halo peaks at 1.2 



X 



10"* erg s"'cm"^sr ' which translates into a dereddened value 
of 4.0x 10"* erg s^'cm^^sr"', ~ 26 times higher than the predic- 
tion. Taking the PSF scattering into account, the 77 = 4 (3) model 
is 4 (14) times below our dereddened upper limit. A similar re- 
lat ion emerges for the upper limit of 1.5 x 10"^ erg s^'cm^^sr"' 
bv lFesen et all (1 19971) . Their result is based on unfiltered spectra 
and must therefore be compared with the dotted lines in Fig. [TT] 
at roughly 0.3 arc min from the edge of the nebula. Including 
reddening, the rj - A {3) model is 4 (10) times below their 
detection limit . For a temperatur e of the shell of 7.5 x 10^ K, 
as assumed bv lFesen et al.l (1 1 9971) . the modelled surface bright- 
ness is ~ 2.7 times higher than for 2 x 10^ K. This brings the 
predicted surface brightness closer to the imaging and spectro- 
scopic upper limits by the same factor. For the 77 = 4 (3) model 
the limits would only be a factor of ~ 1 .5 larger than the emit- 
ted surface brightness. The recombination emission can also be 
enhanced to some degree if the gas is clumpy and has a higher 
gas density. The latter, however, would quickly exceed a sensible 
va lue for the shell mass, m ake the C IV /il550 absorption seen 
bv ISollerma n et all (l2000l) too high, and contradict the findings 
oflSankrit et al.| ( Il998b" 



The slope of the haloes expected from recombination emis- 
sion is close to exponential with 10"°*"* for i] - 3 and 
for T] - A. We note that the latter is coincidentally very simi- 
lar to the slope for the WFI halo caused by PSF scattering. The 
MOSCA data would require steeper density profiles than 77 = 4. 
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R„i„ = 5x10 cm, T, = 2x10 K, M = 4.5 M„ 




0.0 



0.5 1.0 
r [ a r c m i n 1 



Fig. 10. The MOSCA halo profile and the best fit PSF contribu- 
tion (dotted line). The other lines show the Icr uncertainty of the 
fit. 



4.2.2. Dust scattering 

Scattering by interstellar dust can contribute to the shell flux, 
too. The Crab nebula is an isolated object, and its dust den- 
sity is nearly equal to the surrounding ISM, as was found by 
observations from the optical to sub-millimeter \yavele ngths 
(iFesen & Blaidll990l:lGreen et al.ll2004l:lTfemim et al.> 2Q06'). The 
latter authors also showed that small-grain dust is entirely ab- 
sent inside the nebula, and that larger grains amount to only 
10"^ - 10"^ Mq. Hence a normal dust-to-gas rati o can be as- 
sumed and an analysis similar to that in (Osterbrock 19891 chap- 
ter 7) is val id. We adopt an H a luminosity of 7 x 10^"^ erg s"' for 
the nebula dPavidsonll 19871) and take the electron density to be 
3 cm""* close to it, as is appropriate for the t] - 4 model described 
above. In this way one obtains a ratio of emissivities (scattered 
emission on dust to recombination emission) jo/ jua ^ 1- Thus 
dust scattering will contribute about the same amount to the peak 
halo surface brightness as recombination emission. 

The ratio jo/ jua has a radial dependence of <x (7?//?niin)''^- 
For T] - 4- this means that dust scattering will dominate over re- 
combination emission for increasing distances, effectively low- 
ering the steepness of the halo profile. This can be used to dis- 
criminate scattering from recombination emission. Such an at- 
tempt should be complemented by deep observations in [O III], 
as the Crab nebula is bright in this line and dust scattering more 
eflicient for shorter wavelengths. 



4.2.3. Collisional excitation 

One way to significantly enhance Ho- emission is collisional ex- 
citation of Lyp, followed by part of the radiative de-excitation 
going into Ha. This only works for temperatures well above 
10^ K, agreeing with our assumptions of 2 x 10"* K. A neutral 
fraction of hydrogen of 10% is then enough to increase Ha by 
a factor of ~ 10 compared to pure recombination emission. If 
the shell gas was clumpy it is likely that such conditions pre- 




-.5 .5 1 1.5 2 

Distance from edge of Nebula (arcmin) 



1.5 



— 




-1 -.5 .5 1 1.5 2 

Distance from edge of Nebulo (arcmin) 

Fig. 11. Ha halo surface brightness using the model described in 
Sect. 14.21 The dotted lines show the unfiltered total Ha surface 
brightness, whereas the solid and dashed lines show what is cap- 
tured through Ha and continuum ESO/WFI filters, respectively. 
Note how the surface brightness falls off rapidly with distance 
from the edge of the nebula. As discussed in the text, the 77 = 4 
model provides an adequate fit to the shape of the fall-off, but it 
cannot explain the flux amplitude. 



vail. However, our observations and those o f iFesen et alj(ll997l) 
demand collisional excitation to be less dominant, as otherwise 
the halo would have been detected unamiguously. It remains un- 
clear though if this effect will be constant with distance from the 
edge of the nebula, which could alter both the halo's slope and 
amplitude. Detailed modelling is required to test this. 

5. Prospects for a spectroscopic halo detection 

Direct imaging is not the only option for detecting a fast shell. 
iFesen et aiT(ll997h took 9 ksec spectra with a 2.4m telescope of 
regions outside the nebula, but could not detect any emission 
lines. Such emission would only have small Doppler shifts as 
the gas in this area moves mainly perpendicular to the line of 
sight. Another way would be spectroscopy of objects near the 
explosion centre, where the gas moves mainly along the line 
of sight and hence has a very large Doppler shift. In this way 
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Fig. 12. Dereddened spectrum around Ha 10" north of the pul- 
sar In the lower panel we have inflated the flux scale and 
added the rj - 4 model spectrum from Fig. [13] A constant of 
10"'^ erg s^'cm^^sr"' was included, simulating the continuum 
emission from the pulsar wind nebula. Note that the model spec- 
trum has been multiplied by a factor of 100 to be distinguishable 
in this plot. 



ISollerman et al.l (l2000l) detected blue-shifted C IV absorption 
which is currently the best observational evidence for a fast shell. 

Below we evaluate the chances of detecting the Ha footprint 
of a shell towards the centre of the nebula. It would appear as 
emission on both sides of the Ha line, serving as an additional 
verification of an expanding shell. Yet such an observation would 
be made against a complex and bright Ha background, and thus 
the location where the spectra are taken must be picked care- 
fully. A good choice will maximise the contrast as well as the 
spectroscopic separation of the shell and the nebula emission. 

A set of FORS 1 dAppenzeller et al .119981) VLT spectra of the 
pulsar, taken for a different programme, serves as our testbed. 
The total exposure time was 600 seconds through a I'.'O wide slit, 
the dispersion is 2.7 A pixel We identified an area 10" north 
of the pulsar which is largely free of bright filaments (see Fig.|2]). 
The spectrum around the Ha line is shown in Fig. [12] 

We base our analysis on the shell models from Sect. |4.2| and 
their recombination emission. The Ha line profile from such a 
shell depends strongly on the positioning of the spectroscopic 
slit, as gas with different velocities is projected along different 
lines of sight. For example, as the shell has an inner radius Rmin, 
no emission close to the pulsar will be seen for velocities \V\ < 



T \ \ \ I 1 



at the center ij = 4 




-3000 -2000 -1000 1000 2000 3000 



Velocity (km/s) 

Fig. 13. Line profiles for the Ha recombination emission of a 
fast shell, based on the models of Sect. 14.21 Note the gap in halo 
emission between +1680 km s"' for lines of sight close to the 
explosion centre. The smaller the slope rj of the density profile, 
the less pronounced is the intensity fall off with increasing abso- 
lute velocity. The dotted line shows the filter curve of the WFl 
red continuum filter. 

1680 km s"' (+37 A) from the rest wavelength of Ha. This gap 
decreases for lines of sight further away from the centre. At the 
edge of the nebula the halo emission peaks at the rest wavelength 
of Ha. The corresponding spectra are shown in Fig. [13] 

In Fig.[T2]we overplotted the emission expected for the rj - 4 
model over the Ha line observed. The model was increased by 
a factor of 100 compared to the data in order to make it visible. 
Even if collisional excitation was very important and boosted 
the halo's surface brightness by a factor 100, the emission would 
still be undetectable in these spectra. The main problem are the 
neighbouring [N 11] lines at 6548 A and 6583 A corresponding 
to -675 km s"' and +944 km s"', respectively. Their rest wave- 
lengths are close to the expected Ha peak emissions of the shell. 
The motions of the Ha filaments with ~ +1400 km s"' contami- 
nate the halo signal further Increasing the integration times sig- 
nificantly would not help to resolve this issue. We tested other 
areas along the slit as well and they produce no better results 
than this one. We therefore conclude that a search for Ha shell 
emission inside the Crab nebula is unfeasible. 



6. Discussion 

6. 1 . PSF scattering and haloes around galaxies 

Our findings concerning the PSF scattering are applicable to 
other situations where very faint haloes are searched for around 
bright extended targets. The effect is easily overlooked, as the 
PSF is often thought of in scales of arcseconds rather than of ar- 
cminutes. The PSF can be very efficient in scattering light from 
far inside a bright target into its outskirts. 

Observations of the extraplanar gas around edge-on galaxies 
are particularily affected by this, as galaxies can have high cen- 
tral surface brightnesses and small minor axes. Wuetal. (200l) 
recognised this effect in their work on NGC 4565 and con- 
structed large PSF models based on separate observations of iso- 
lated bright stars. Their models extend to 28' (corresponding to 
1000 pixels), and they fin d that PSF scatte ring in their data is 
measurable but negligible. IWu et all (|2002l) are cited by several 
authors working on similar subjects, yet PSF scattering appears 
to be largely ignored. 
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IZibetti et alj ( |2004|) take this effect into account in their anal- 
ysis of the stacked halo of 1047 edge-on galaxies taken from 
SDSS. They find a maximum contamination of about 30%. 
However, their PSF extends only to about 20 pixels (8") and 
is then extrapolated by an exponential function. Ide Jond (|2008|) 
investigated this in detail and reconstructed the SDSS PSF out 
to 180 pixels. Beyond that, the PSF is extrapolated with a power 
law index of -2.6 as obtained from the inner 40 to 18 pixels. 
Using a full 2D convolution of the stacked galaxy image. lde Jond 
(l2008h shows that the PSF contributes 5 0% in g and r , and u p to 
80% in /-band to the signal measured bv lZibetti et"an(l2004 . 

Even though this is an entirely different field of astrophysics, 
it suffers from the same observational effect. In our case the sit- 
uation is rather extreme as the Crab nebula is of effectively con- 
stant surface brightness over its 5' x 7' area when smoothing out 
the filaments. Thus PSF scattering is ver y effective as t he en- 
tire object contributes. The conclusions of Ide Jon^ (|2008|) are in 
excellent agreement with our own: 

- Extrapolating a PSF to large radii can be very misleading as 
a small change in the wing amplitude can have a large effect 
on the amount of scattered light. 

- The PSF wings should ideally be constructed from extremely 
bright isolated stars and be joined with a core determined 
from bright but unsaturated stars. 

- The amount of scattering is very difficult to judge based just 
upon the PSF profile. Only a full 2D convolution will reveal 
its overall effect. 

6.2. Prospects for future observations 

Our images are strongly affected by PSF contamination. This 
also applies for similar attempts at other telescopes, and can only 
be overcome by large and accurate PSF models. These must be 
created by means of a very bright star that by far dominates all 
other flux in the image. Ideally, the star should be positioned on 
the same detector area as the target. It should be bright enough 
so that the PSF can be modelled without using azimuthal averag- 
ing, which destroys asymmetries in the PSF wings. Such a star 
should ideally be observed close to the target and in the same 
night, bracketing the target observations in time. The detector 
must be given sufficient time for recovery after such observa- 
tions so that charge persistencies can decay. Dark conditions are 
required to keep the sky background as low as possible. 

Observations through a red continuum filter would signifi- 
cantly minimise the contamination by scattered light. We inves- 
tigated the efficiency of this approach based on the WFI con- 
tinuum filter As can be seen from Figs. [TTI and [T3l the flux ex- 
pected in that filter would be 5 (8) times lower for the 77 = 3 (4) 
models as compared to Ha. This will be very difficult to record 
on 8m class telescopes. On the other hand, the combined emis- 
sion from the processes discussed in Sect. 14.21 or from its con- 
stituents alone, would become visible in several hours Ho- expo- 
sure with e.g. FORS2 on the VLT, provided that good PSF mod- 
els are available. For an optimum comparison of such data with 
model predictions one would subtract the 2D PSF-convolved im- 
age from the original image, and then integrate over the halo in 
concentric rings around the explosion centre. 

On the spectroscopic side, the most promising approach is 
to search for Doppler shifted absorption features in a relatively 
uncontaminated spectral region against the central part of the 
nebula. However, care has to be taken in the interpretation of 
the result, as a chimney-like structure extending towards us may 
contain gas at velocities similar to those in the shell. This can be 



tested by searching for fast (~ 2000 km s ') blueshifted [O III] 
A4959 emission as [O III] is bright in the Crab chimney. This 
would lead to emission at around 4925 - 4935 A and could be 
difficult to disentangle from He 1 492 1 . Adding information from 
other emission lines seen in the chimney may be useful. 

Finally, deep 21 cm radio observations outside the nebula 
may be useful. The main difficulty here is the uncertain degree 
of ionisation of hydrogen in the shell. If collisional excitation of 
H-alpha is important, then a good fraction of neutral hydrogen 
is implied which may show up in 21 cm line emission. Hence 
observations of the Ha surface brightness of the shell can help 
to estimate the radio flux. 

7. Summary 

We have presented WFI and MOSCA images in Ha to search 
for a high velocity shell around the Crab nebula. Our detec- 
tion limit in both data sets is 5 x 10"** erg s~'cm ~^sr~', three 
times deeper than the previously deepest study by iFesen et all 
(Il997h . However, our measurement of an Ha halo is corrupted 
by PSF contamination which can account for all of the signal. 
A real halo with an observed (i.e. reddened) peak brightness 
of 2 X 10"^ erg s^'cm^^sr"' could still be accomodated within 
the error bars. Thu s we could not improve on the limits set by 
iFesen et al. 1 (Il997h . Deep imaging with 8m class telescopes, to- 
gether with good PSF models, could yield progress in this field. 

We discuss three different processes that could power Ha 
emission from a fast shell: recombination emission, dust scatter- 
ing, and collisional excitation. All but the latter yield flux levels 
about one order of magnitude below the one observed. As the 
halo flux in two independent imaging data sets is well described 
by PSF scattering, we can rule out collisional excitation to play 
a dominant role in the shell gas. 

We argue that a spectroscopic detection of Ha emission to- 
wards the centre of the nebula is unfeasible. The emission is very 
weak and buried in the complex spectral environment of the Ha 
and [N II] lines. The prospects of detecting blueshifted absorp- 
tion features in the pulsar or other background s ources are signif- 
icantly better and were already demonstrated bv lSoUerman et all 
( l2000h . A deeper spectrum outside the Crab nebula could still be 
very useful, as the recombination level is below the upper value 
from iFesen et all (Il997h . However, the amount of PSF scatter- 
ing needs to be quantified also for spectroscopy. This problem 
could be avoided by looking only at emission further away than 
about 2200 km s ' from Ha. However, the flux expected for den- 
sity profiles steeper than ?/ ~ 2 is so low that a positive detection 
appears unfeasible. 
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Appendix A: MOSCA peculiarities 

We discovered that the gains of the individual chips in the 
MOSCA array depend on the illumination level. The gain ra- 
tios are different by up to 20% in flat fields with approximately 
30 kADU as compared to the observations with very low back- 
ground level (~ 100 ADU). The same behaviour was observed 
in data taken through two other filters about two months later. 

The THELl pipeline assumes that all chips of a mosaic cam- 
era are brought to the same gain during flatfielding. This is not 
the case for MOSCA. We circumvented the problem by switch- 
ing ofi" THELI's gain adjustment and processed and stacked the 
data from each CCD independently. The 4 stacked images were 
then run again through the astrometric and photometric part of 
THELl, being interpreted as 4 dithered exposures coming from 
a single-chip camera. Given the overlap sources, relative photo- 
metric zeropoints could be calculated and thus the gain differ- 
ences were removed. Upon comparing the fluxes of 74 common 
stars in the MOSCA and WFl data, we find a scattering of 0.08 
mag without systematic trends as a function of position in the 
MOSCA mosaic. 



